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We have shown, that the possibihty of the existence of the mixed phase of the non-uniform chiral 
(NCh) and the color superconducting (2SC) ground state depends significantly on the choice of the 
parameters and type of the regularization scheme. Our calculations indicates, that in the 3d cut-off 
regularization scheme, the mixed region of the NCh and the 2SC phases exists for a broad set of 
NJL model parameters. However, in the Schwinger regularization scheme, if parameters are set to 
the vacuum values of /^, m-^ and < qq >, then, the mixed region of the NCh and the 2SC phases 
does not exists. 
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Q^. I. INTRODUCTION 

^ ,: 

C/j . The goal of this study is to examine the influence of the current quark mass on the phase structure of the strongly 

\^ I interacting quark matter. The part of the QCD phase diagram that was attracting peoples attention for many years 

is a region of relatively low temperatures and moderate baryon densities. In this region, a variety of phases are 

possible. Especially, a competition between the chiral condensate and the diquark condensate takes place [l| and 

f\, both these condensates might coexist. Moreover, since a spatially non-uniform chiral ground state was suggested 

I ' [3] J the situation is even more complicated. Many studies was devoted to this issue and most of them was done in 

Mh| the chiral limit. Taking into account both a non-zero current quark mass and a chiral density wave ansatz, makes 

~^ ' calculation difficult. The interesting method to deal with that mathematical inconvenience was suggested recently 

[si |j|. The clue is to expand the quark propagator in a perturbative series. We will follow that method, working 

with the Nambu - Jona-Lasinio model, with two flavors and three colors. It is now well known, that within the 

Cn ' mean field approximation in the chiral limit, the chiral density wave ansatz guarantees a lower minimum (in the 

moderate baryon density region) than a uniform chiral phase. But what happens if we include the non-zero current 

quark mass? In Ref. [3] it was shown, that in the Schwinger regularization at zero temperature, introduction of the 

\^^ non-zero current quark mass does not exclude the possibility of existence of the spatially inhomogeneous ground 

pO ' state. However, there are open questions. First of all, whether, away from the chiral limit, the non-uniform chiral 

phase manifests also in different types of regularizations. This is a crucial aspect, especially when doing model 

dependent calculations [5|. Secondly, there is a question, how the non-zero temperature does influence the phase 

p^ ' diagram? Finally, we would like to resolve whether there is a possibility of the existence of the mixed region of the 

2SC color superconducting and the non-uniform chiral phases. 
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U{ ■ II. MODEL 

C^ , 

The starting point is based on the Nambu - Jona-Lasinio model with two flavours and three colors. Because we 
will often refer to Ref. [3|, we must notice that in the present case, the initial lagrangian density differs by adding a 
charge conjugate fields tpc and by including a diquark terms. 



H = [^{i^di, + A<7o - m)'(/' -t- TpciiYd^, + ^70 - m)'ijjc 

J X 

+ G [(ViV)' + i^ll5T^f] + G'{^cll5T2X^i'){^il5T2X^^c)i (1) 

where ip is the quark field, tpc — C'lp^ is the conjugate field and jjl is the quark chemical potential. We assume that 
current quark masses are equal, ruu = rrid = m. The value of the mass m is fitted with the other parameters of the 
model. The color, flavor and spinor indices are suppressed. The vector r is the isospin vector of Pauli matrices and 
A"^, A — 2,5,7 are three color antisymmetric SU{3) group generators. The integration J^ = J^ dr J d^x, where /3 

is the inverse temperature and the derivative operator is d^ — {idr, V). Two coupling constants G and G" describe 
interactions which are responsible for the creation of quark- antiquark and quark-quark condensates respectively. Both 
couplings are treated as independent. There is also an additional parameter A which defines the energy scale below 
which the effective theory applies. It is introduced through the regularization procedure. We are working in the mean 



field approximation within the ansatz 
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which describes three possible phases: the chiral uniform phase Ch {q — 0, M ^ 0, A — 0), the non- uniform chiral 
phase NCh (<f 7^ 0, M ^ 0, A = 0) and the color superconducting phase 2SC (g = 0, Af = 0, A 7^ 0). Introduction of 
the non-zero current quark mass results in a fact that a chiral phase order parameter is always larger than the mass 
m, and the chiral symmetry is only partially restored. The presence of the non-zero current mass, causes also, that 
the 2SC phase always coexists with a small quark condensate. We will perform the calculations in a Nambu-Gorkov 
basis X"^ = (■0,'0c) \3i- To deal with the dependence of the quark propagator on the space coordinates, we make a 
standard chiral rotation. The mean-field partition function has a form 
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and 



To determine the ground state of the system we calculate the thermodynamic potential 

Expanding the termodynamic potential up to the terms of order Vm we decompose f2 into the two parts [3| 
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Using the standard method one can calculate the thermodynamic potential ilo |8| 



0{V^. 



(9) 



no 



(Mt - mf 
AG 

(2^ 

d3fc 



^-Sl.(0^--^-S/(0(-S-^) 



y 2T I -— 



In 1 -f exp 



In I 1 -I- exp I — 



^/{ElTW+W\ 



T 



T 



In 1 + exp — 



s=± 

In 1 + exp 

Es — /i 



T 



T 



where 
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The term SH, gives a non-zero contribution to the potential only if g ^ 0, and in a momentum-frequency basis, has a 
form [3| 
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The result of Matsubara sum is rather complicated and more detailed study can be found in the Appendix 
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III. REGULARIZATION PROCEDURE 



The Nambu - Jona-Lasinio model is a nonrenormalizable theory. That is why, the regularization procedure is needed 
to clearly define the model. We use the two types of regularization procedure: 3d cutoff and the Schwinger (proper 
time) scheme. The temperature dependent part of the J7o is finite, but the zero temperature limit of the potential D,q 
contains divergent terms. We rewrite ^o(T = 0) in a form that is more convenient to regularize 
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Only the last two terms of equation (|T4)) are divergent. The infinite contribution to the SQ comes from the following 
terms 
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3d cutoff regularization 

We regularize divergent integrals by introducing a cutoff A, k^ < A. 



• Schwinger (proper time) regularization 

This method is more subtle Q. At first, we decompose divergent integrals into the terms that depends on 
temperature (these terms are finite) and the remaining part. The finite terms contains factors of a type (1 + 
exp{~l3E^ gj)^^ . With momentum tending to the infinity, above mentioned factors tends to 1. In the example, 
with this method, integral ([T7| has a form 
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The last term of equation (jlSp is finally regularized in the Schwinger scheme 
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3d 



A 


0.653 1.086 


GA^ 


2.1 3.68 


Mo 


0.313 0.2 


m 


0.005 0.0049 



{-<uu>)^ 0.25 0.25 



TABLE I: Numerical values of the regularization parameters in GeV (type 1) 



3d 



A 


0.619 0.635 


GA^ 


2.24 7.53 


Mo 


0.35 0.38 


m 


0.0053 0.0153 



{-<uu>)3 0.245 0.17 



TABLE IL Numerical values of the regularization parameters in GeV (type 2) 



IV. RESULTS 



The parameters of the NJL model are G, G", A and in. As we shall see, the particular choice of the above mentioned 
parameters can have a decisive influence on the conclusions following from the NJL model. The physical quantities 
in the vacuum that are used to fix G,A and m are: /^r = 93 MeV, m.^. = 135 MeV and < uu >=< dd >= -(250 
± 50 MeV) 3 [10|. Because of the wide experimental range for the value of the quark condensate density, one can 
neglect < qq > and then, A becomes a free parameter yy|. To make the discussion clear, we define the parameters 
of type 1 (when the < qq > value is taken into consideration) and of type 2 (when A is a free parameter) . The fitted 
parameters are collected in the Tables I, II, where also the value of the constituent quark mass in the vacuum Mq 
is given. Coupling constant G' can not be related to any known physical quantity. The larger the value of G' the 
lower the value of the fio in the minimum corresponding to the 2SC superconductor. Because we want to resolve 
whether the mixed region of the NCh/2SC phase occurs when to 7^ 0, we set G' = 0.5 G. That is a lower bound 
for the coupling constant G' considered among the vast literature on the color superconductivity [a, [I4I . With the 
fixed parameters G, G', A and to, the thermodynamical potential 51 is a function of Mt,q, |A|, chemical potential n 
and temperature T. We begin ours analysis with a zero temperature limit. For a constant /i and T, we numerically 
determine a global minimum of J7 with respect to the free parameters. We restrict our analysis to the region of /z, 
where numerically obtained Mt is greater than mass m. 



A. 3d cutoff regularization 
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FIG. 1: The values of Aft, q and A in the 3d cut-off regularization scheme for T=0 (GeV units). The global minimum related 
to the mixed phase NGh/2SC is plotted, (type 1 parameters) 
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FIG. 2: The phase diagram in the /i - T plane, (type 1 parameters) 



Type 1 parameters 



In the Fig. 1, dependences of Mt^q and |A| on /i at T = are shown. The spatially inhomogeneous quark phase 
does exist. There is a global minimum corresponding to the mixed phase of non-uniform chiral (NCh) and color 
superconducting (2SC) phases. At zero temperature, its value is slightly lower than the value of the local minimum 
corresponding to the mixed phse of 2SC and uniform chiral (Ch) phases, in the interval of /z from 0.322 GeV to 0.420 
GeV (Fig. 1). As we can expected [8|, non-zero temperature does not favours a non-uniform condensate. The phase 
diagram in the /i - T plane is shown in Fig. 2. For temperatures between and 20 MeV, the interval of existence of the 
inhomogeneous phase is practically constant, but the region of a coexistence of NCh and 2SC phases is getting smaller 
with growing temperature (Fig. 2, 3). Above T = 16 MeV, spatially inhomogeneous quark condensate exists without 
a diquark condensate, and region in which the global minimum is associated to NCh is systematically decreasing. 
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FIG. 3: The values of Mt, q and A in the 3d cut-off regularization scheme for T=0.011 (GeV units, type 1 parameters). 



2. Type 2 parameters 
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FIG. 4: The values of Mt, q and A in the 3d cut-off regularization scheme for T=0 (GeV units). The global minimum related 
to the mixed phase NCh/2SC is plotted, (type 2 parameters) 



The wide experimental range for the quark condensate density allows us to shift the value of A. Moving A above 
the previous value 653 MeV results in a decrease of a value of the coupling constant G. For the lower values of G, 
minimum associated to the NCh phase is getting more shallow, that case will also be discussed later. Hence, we will 
rather reduce the value of A. On the other hand, A is an upper boundary for the energy scale, so the possible interval 
for the values of A is narrow. We set GA^ = 2.44, A = 619 MeV, Mq = 350 MeV and m = 5.3 MeV. With this choice 
of parameters, (— ^ < qq >)3 — 245 MeV and it is still within the experimental range for the value of the quark 
condensate density. Fig. 4 presents the dependences of Mt,q and |A| on /i at T = 0. The spatially inhomogeneous 
quark phase exists in a region of /i from 340 to 440 MeV. In that region, the 2SC color superconducting phase coexists 
with the NCh phase. As one can expected, with growing T, above mentioned region is getting narrower and the value 
of the color superconducting order parameter |A| is decreasing (Fig. 5). Above T = 20 MeV, there is no mixed region 
of the NCh/2SC phase and these phases exists separately. 
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FIG. 5: The values of Mt, q and A in the 3d cut-off regularization scheme for T=0.016 (GeV units). The global minimum 
related to the non-uniform chiral solution is plotted. 



B. Schwinger regularization 
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FIG. 6: The values of Mt and A in the Schwinger regularization scheme for T=0 (GeV units). 



Type 1 parameters 



In the proper time scheme at T = 0, the mixed phase of the NCh and 2SC phases is only a local minimum, and 
the values of the vector q corresponding to that local minimum are unphysical {q is much greater than the A scale). 
In Fig. 6, there are presented dependences of Mt and |A| on /i at T = 0. Even at T = 0, we observe a smooth 
dependence of Mt on a chemical potential. Mt decreases from its vacuum value, across the narrow range of the 
rapid variability, down to the 5 MeV. In this regularization scheme, the values of |A| are of order of the magnitude 
smaller than in the cutoff scheme. That fact was already observed Q. For the values of T above 5 MeV, the diquark 
condensate vanishes. 



2. Type 2 parameters 



The limit values of the NJL parameters for which (-5 < qq >) 3 > 200 MeV are GA^ = 4.58, A = 760 MeV and 
m = 12.6 MeV. However, at these limit values, there is no non-uniform phase. Guided by the findings of [3|, we set 
Gh? = 7.53, A == 635 MeV, Mq = 380 MeV and m = 15.3 MeV. With this choice of parameters, (-i < qq >)^ = 170 
MeV and it is already beyond the experimental range for the value of the quark condensate. As can be seen in Fig. 
7, the mixed region of the non-uniform chiral/2SC superconductor ground state appears from /i — 355 MeV. Above 
the value of quark chemical potential equal 450 MeV, the value of the wave vector q is greater than A. Therefore, 
above these limit, discussed model is unphysical. 
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FIG. 7: The values of Mt, q and A in the Schwinger regularization scheme for T=0 (GeV units). GA^=7.53, A=0.635 GeV, 
m=0.0153 GeV 



C. Discussion 

To compare the behavior of the potential 57 in different regularizations and at the different choice of parameters, 
we take the factor A^ ahead of the Vl and then we work with the dimensionless quantities. The 50, contribution to 
the thermodynamic potential is proportional to the ratio m/ K. When the value of q is non-zero, 5^1 is positive and 
therefore, the value of Q, is higher than the value of JIq. As a consequence, 5Vt term does not favour the possibility of 
the existence of the NCh phase. One might ask a question, why so small mass m (in comparison with other energy 
scales) has so significant infiuence on the J7. In the 3d cutoff scheme with the parameters of type 1, the ratio 7n/A 
= 0.0076 and with type 2, m/ K = 0.0085. Similarly, in the proper time scheme, with the parameters of type 1, the 
ratio m/K — 0.0045 and with type 2, m/A — 0.024. It appears, that the difference between the value of the fio in the 
global minimum corresponding to the NCh/2SC phase and the value of the Oq in the local minimum corresponding 
to a uniform chiral phase is also of order of a few percent. Therefore, even so small mass m can moves the global 
minimum. There is also another question, why the inhomogeneous ground state solution appears for the rather higher 
current masses. Of course, the parameters of the NJL model are linked with each other. It turns, that a larger m 
requires a larger coupling constant G. With increasing G, the distance between the value of the fig in the global 
minimum corresponding to the NGh/2SG phase and the value of the f^o in the local minimum corresponding to a 
uniform chiral phase is more distinctive. In Fig. 8. we show the behavior of the values of the full potential fi in the 
minima corresponding to the NCh/2SC and Ch/2SC solutions as a function of mass m in the 3d cutoff scheme. The 
upper limit on the value of mass m is m k 5.45 MeV when the vacuum value of Mq exceeds 400 MeV. On the other 
hand, below the value of mass m k, 4.7 MeV the non-uniform solution disappears. 



CONCLUSIONS 



We have shown, that the possibility of the existence of the mixed phase of the non-uniform chiral (NCh) and 
color superconducting (2SC) ground state depends significantly on the choice of the parameters and type of the 




m [MeV] 

FIG. 8: The values of the Q in the mmima corresponding to the NCh/2SC phase and to the Ch/2SC phase (dashed line), as a 
function of mass m in n — 360 MeV (3d cut-off, type 1 parameters) 

regularization scheme. Our calculations indicates, that in the 3d cut-ofF regularization scheme, the mixed region of 
the NCh and 2SC phases exists for a broad set of NJL model parameters. However, in the Schwinger regularization 
scheme, if parameters are set to the vacuum values of /tt, Wtt and < qq >, then, the mixed region of the NCh and 2SC 
phases does not exists. In the moderate baryon density region, there is a local minimum associated to the NCh and 
2SC phases, but the ground state is realized by a chiral uniform, or in a small range of temperatures above zero by 
the uniform chiral (Ch) and superconductor (2SC) phases. If we release the constrain on the parameters coming from 
the value oi < qq >, then, the NCh solution might be an absolute minimum of the thermodynamical potential in the 
moderate baryon density region. The result, that is irrespective for the parameters and regularization method choice 
is that, if for some chemical potential (at T = 0) there is a minimum associated to the NCh/2SC phase, its value is 
always lower than the value of the minimum associated to the spatially inhomogeneous quark condensate with |A| 
equal to zero. Because we chose G'= 0.5 G, that is a lower bound for the estimated value of a coupling G", it can be 
assumed, that if the NCh condensate exists at zero temperature, it exists together with a diquark condensate. With 
growing temperature, a diquark condensate starts up at a greater n and the region of a coexistence of the NCh and 
2SC phases disappears. 



VI. APPENDIX 



In the expression ([121 
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the trace is taken over the Nambu-Gorkov, color and flavor indices. The sum is taken over the Matsubara frequencies. 
The diagonal elements of the Sq operator Eq. ([5]) in the momentum frequency bassis have a form 
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where indices a,^ refers respectively to the flavor and Nambu-Gorkov space and are equal 1, 2. Index /3 refers to the 
color space and takes the values 1, 2, 3. We calculate that trace in the eigenbasis of a matrix 70(^7 — Mt ± §7573). 
The eigenvalues of a matrix joikj — Mt ± §7573) to the eigenvectors tpi are Xi ~ ±E± (fTT|) . For example, the trace 
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of a one of these diagonal elemsnts is 
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Finally we get 
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Expression p^ . is obtained by the contour integration of expression ([^^ JJ and by a summation over the whole set 
of indices. 
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